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Abstract. We present a new characterization of partially co- 
herent electric and magnetic wave vector fields. This char- 
acterization is based on the 36 auto/cross correlations of the 
3 + 3 complex Cartesian components of the electric and mag- 
netic wave fields and is particularly suited for analyzing elec- 
tromagnetic wave data on board spacecraft. Data from space- 
craft based electromagnetic wave instruments are usually 
processed as data arrays. These data arrays however do not 
have a physical interpretation in themselves; they are simply 
a convenient storage format. In contrast, the characterization 
proposed here contains exactly the same information but are 
in the form of manifestly covariant space-time tensors. We 
call this data format the Canonical Electromagnetic Observ- 
ables (CEO) since they correspond to unique physical ob- 
servables. Some of them are already known, such as energy 
density, Poynting flux, stress tensor, etc, while others should 
be relevant in future space research. As an example we use 
this formalism to analyze data from a chorus emission in the 
mid-latitude magnetosphere, as recorded by the STAFF-SA 
instrument on board the Cluster-II spacecraft. 



1 Introduction 



When analyzing time varying electric and magnetic vector 
field data from spacecraft, it is common to construct a 6 x 6 
matrix from the complex vectors (E,cB), sometimes jointly 
called a sixtor. In a Cartesian coordinate system, this energy 



density matrix (in SI units) can be written 
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This electromagnetic (EM) sixtor matrix has in vari- 
ous guises, such as wave -distribution functions (WDF) 
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dStorev and Lefeuvrd 119741) and so on, been useful in the 
analysis of EM vector field data from spacecraft, for instance 
on the Cluster and Polar missions. The second order co- 
herency matrix is important from a statistical viewpoint since 
it completely describes a wide sense stationary vector signal. 
From a physical point of view it is important since energy 
density, EM wave polarization, and similar quantities can be 
derived from its components. 

The EM sixtor matrix can be seen as a generalization of 
the coherency matrix in optics, which is usually a 2 x 2 Her- 
mitian matrix, describing the transverse field. The coherency 
matrix description is convenient as a data storage format but 
in practice, it is more common to instead use the four Stokes 
parameters, as they are physically more intuitive. 

One may ask what the Stokes description for the full EM 
field would be? One way is to decompose the 6x6 coherency 
matrix in terms of a complete b asis set of unitary matrices, 
as per lSamson and Olsenl Jl98dl) . The problem with this ap- 
proach is that such a decomposition is not unique; there is an 
infinite number of unitary bases. 

In this paper we introduce a unique set of parameters, anal- 
ogous to the Stokes parameters, but generalized to the full 
electric and magnetic wave fields. We call these parameters 
the Canonical Electromagnetic Observables (CEO), due to 
their uniqueness, which comes from the fact that they are ir- 
reducible under Lorentz transformations. Furthermore, they 
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are not merely "parameters", they are proper space-time ten- 
sors. 

Some examples of the CEO in vacuum are the EM en- 
ergy density £o(|E| 2 + |B| 2 )/2, and the energy flux density 
(proportional to the Poynting vector) Re [E x B*] /Zq, where 
Zq = y/no/eo is the vacuum impedance. These quantities 
can be quite easily identified from the sixtor matrix by in- 
spection, for instance its trace is the energy density. Other 
CEO are not so easily identified. 

The procedure to calculate the fundamental CEO, in 
terms of six real irreducible space-time 4-tensors, as well as 
the corresponding three-dimensional representation will be 
given in Section 2. A two-dimensional CEO representation 
is given in Section 3. 

1 . 1 Covariance in space physics 

Perhaps even more important for space borne observations: 
the EM sixtor matrix is not covariant according to the re- 
quirements of special relativity. Spacecraft are constantly 
moving and often spinning observation platforms. EM wave 
measurements become Doppler shifted and data must often 
be "despun". For the four 3x3 sub-matrices: £oE®E*, 
E(g)B*/Zo, B(g)E*/Z , and B®B*/Mo, where ® denote the 
outer (tensor) product, despinning is straight forward by ap- 
plying rotation matrices R from left and R T from right, e.g. 
E'(g)E' = RE (g)ER T . To rotate the full EM sixtor matrix, 
similar operations must be performed four times. This is 
awkward and the resulting 6x6 matrix is still not covariant. 
For EM wave measurements in space plasma the last remark 
can be crucial. 

A Lorentz boost is the translation from one Lorentz frame 
to another one moving at velocity v. A Lorentz boost does 
not necessarily imply relativistic speeds, which is a common 
misconception; and therefore it do not by itself preclude what 
is typically associated with relativistic effects. It is simply a 
quite general recipe to make two different observers agree 
on a physical observation. The Lorentz boost of the EM field 
vectors can be written E' = y(E + v x B) and B' = y(B - 
v x E/c 2 ), where y = 1/ a/1 — v 2 /c 2 . As a matter of fact, 
the Lorentz boost is the essence of the well-known frozen- 
in field line theorem^ from magnetohydrodynamics (MHD); 
a theory which is commonly used to model the solar wind 
plasma. In a plasma, relativity comes into play at very a 
fundamental level since the electromagnetic (Lorentz) force 
dominates the vast majority of all plasma interactions. 

Another example illustrates the problem to separate time 
(frequency) and space (wave vector) in EM wave observa- 
tions on board a spacecraft. Assume that we observe a wave 
mode which is described by an angular frequency co and 
wave vector k. We can write this as a 4-vector (o,ck). Let's 



make a Lorentz boost in the v direction: 



'ifE + vxB^Oina plasma, the magnetic field lines change 
as though they are convected with velocity v, i.e., they are frozen 
to the plasma flow. This is the frozen-in field line theorem of ideal 
MHD. 



Co' = y((B k v 
ck' = ck 



V( ck - v ) 



yco 



(1) 

(2) 



What happens now for a stationary (DC) field structure mov- 
ing with the solar wind plasma? We then have 0=0 and 
|ck| ^ 0. For a satellite moving with velocity v relative to 
the DC field structure, it is justified to set y « 1 (the solar 
wind speed seldom reaches more than 900 km/s and using 
this value we obtain y w 1.0000045 > 1); Eqs. © and © 
are then reduced to 



CO 
ck' 



k v 

ck 



(3) 
(4) 



We can see that the DC field structure is not Lorentz con- 
tracted appreciably at this low velocity, k' w k. However, 
there is a dramatic change in the observed frequency, which 
for a head-on encounter with the structure is registered as 
co' ~ kv rather than zero. The observed frequency is propor- 
tional to the dimension of the structure, which we take to be 
in the order of one wavelength, X —2n/k. Taking v = 900 
km/s a 900 km DC field structure would now register as 1 
Hz, a 90 km structure as 10 Hz, and a 9 km structure as 100 
Hz, etc. 

These simple examples clearly show that a space-time (co- 
variant) description is necessary even if y > 1 . The fre- 
quency (time) and the associated wave vector (space) can not 
be treated separately but must be considered together, as a 
space-time 4-tensor. 

The Maxwell equations are inherently relativistic and can 
easily be put into a covariant form using 4-tensors. From 
a theoretical point of view, this fact alone provides a very 
good argument why one should try to express also the sec- 
ond order properties of the EM fields using a covariant for- 
malism. This was recentl y carried out by the authors and 
published in a recent paper lCarozzi and Bergman d2006l) . In 
this paper we introduced a complete set of space-time ten- 
sors, which can fully describe the second order properties of 
EM waves. We call this set of tensors the Canonical Electro- 
magnetic Observables (CEO); in analog y with Wolf's analy- 
sis of the Stokes parameters IWolddl954h . We suggest that the 
CEO could be used as an alternative to the EM sixtor matrix. 
Not only are the CEO covariant, but they are all real val- 
ued and provide a useful decomposition of the sixtor matrix 
into convenient physical quantities, especially in the three- 
dimensional (3D), so-called scalar-vector-tensor (S VT) clas- 
sification; see section [2721 The CEO have all dimension en- 
ergy density but have various physical interpretations as will 
be discussed in what follows. 
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CEO 


Rank 


Proper + 


Number of 


4-tensor 


(Symmetry) 


Pseudo — 


observables 


c+ 





+ 


1 


C- 







1 


T 


2(S) 


+ 


9 


U 


2(S) 




9 


Q 


2(A) 




6 


W 


4(M) 


+ 


10 



Table 1. CEO in space-time classification, i.e., 4-tensor notation: 
1 + 1+ 9 + 9 + 6+10 = 36 observables. 



2 Canonical Electromagnetic Observables 

The CEO set was derived from the complex Maxwell field 
strength F a $ . Other possibilities, su ch as using the 4- 
potential A^ 1 or using a spinor formalism lBarul ( 1980h . were 
considered but discarded due to their lack of physical con- 
tent. The 4-potential is not directly measurable and it is 
furthermore gauge de pendent. Spinor formalism has been 
proved possible to use Sundkvisl (2006) but we believe the 
space-time tensor formalism to be more intuitive and conve- 
nient to use. 

In the quantum theory of light, observables of an EM field 
are ultim a tely c onstructed from a complex field strength; 
see 



Worn (119541) . The simplest of these observables are 
sesquilinear-quadratic (Hermitian quadratic) in F a P, i.e., 
they are functions of the components of F a ^F^ s , which is 
a 4-tensor of rank four. Here we have chosen to denote the 
complex conjugate of the field strength with a bar over the 
field symbol in order not to confuse it with the dual field 
strength, which we denote by a superscript star to the left of 
the field symbol. We showed that it was possible to decom- 
pose F a PpY s into a unique set of tensors, the CEO, which 
are real irreducible under the full Lorentz group. We shall 
not repeat the derivation here but will instead discuss the 
space-time (4-tensor) and three-dimensional (3-tensor) rep- 
resentations of the CEO. 

2.1 Fundamental space-time representation 

In terms of the Maxwell field strength F a P, the 
CEO are organized in the six real irreducible 4-tensors 
C+,C-,Q a Pj a P,U a P, and W a ^ 5 . This is the fundamen- 
tal space-time representation of the CEO; their properties are 
listed in Table 12. ll 

The CEO 4-tensors are defined as follows: the two scalars 
are the vacuum proper- and pseudo-Lagrangians, 



C :={F ap F a P-*F ap *F^)/2, 



C-:=(F a p*F a P+*F a pF a t) f 2. 



(5) 

) /-■ (6) 
respectively, where we have used the dual of F a P defined as 

* F ap^y a pr6 FrS= y^ sF rs_ (7) 



The three second rank tensors consist of the two symmet- 
ric tensors 



(8) 
(9) 



and the antisymmetric tensor 

gap . =i fpa^pMS _ *py F ^ -2C + -n a P) 12. 



(10) 



The symmetric second rank tensor T a P is the well-known 
EM energy-stress tensor, which contains the total energy, flux 
(Poynting vector), and stress (Maxwell stress tensor) densi- 
ties. The other two second rank tensors, U a $ and Q a P, re- 
spectively, are less well-known. The symmetric U a P tensor 
is similar to T a P in that it contains active energy densities but 
in U a & these densities are weighted and depend on the the 
handedness (spin, helicity, polarization, chirality) of the EM 
field. Therefore, we have chosen to call them "handed" en- 
ergy densities. The anti-symmetric tensor Q a P on the other 
hand is very different in that it only contains reactive energy 
densities, which are both total (imaginary part of the complex 
Poynting vector) and handed. 

The fourth rank tensor is 

w aPy8 ._ fpaPpyS _ *pal3* F y8\ ij _ 2i g[a[S r j y]P] 



(ID 



^a[S^y}P _ ^ c _ £ aPyS 



where the square brackets denotes 
tion over the enclosed indices, e.g 



antisymmetriza- 

T^grlP = 

ifpccdgYP _p<xygBp\ an( j nes t e( j brackets are not 
operated on by enclosing brackets, e.g., T^ s g"^^ = 
\ (T a5 grt - T a Yg 5 P - 7^ V" + T^g Sa ) . It fulfills the 

symmetries W a ^ s = W^ 5 = W a ^ = W^ Sa P and 
W a[PyS] _ o 

This real irreducible rank four te nsor, Eg. (fTTTl, was 



discov ered by u^| and published in Carozzi an d Bergman 
d2006h . and is still under investigation; it is an extremely in- 
teresting geometrical object, having a str ucture identical t o 
the Weyl tensor in general relativity; see IWeinberg ( 1972 ). 
We have found that it contains a four-dimensional general- 
ization of the Stokes parameters, as will be demonstrated in 
section[3]for the two-dimensional (2D) case. It contains both 
reactive total and reactive handed energy densities. 

2.2 Three-dimensional representation 

The fundamental space-time 4-tensor CEO can be written 
in terms of the three-dimensional E and B vectors, i.e., 3- 
tensors. This is convenient because it allows us to use in- 
tuitive physical quantities. To systematize the 3D repre- 
sentation of the CEO, we will use a physical classification 



2 To the best of our knowledge, the W a ^^ tensor has never be- 
fore been published in the literature. 
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Scalars 


Vectors 


Tensors 


Active Total (AT) 


u {S AT } 


P {V A r} 


T {T Ar } 


Active Handed (AH) 


w {S AH } 


V {Yah} 


U {T AH } 


Reactive Total (RT) 


1 {Srt} 


R {\ RT } 


X {T RT } 


Reactive Handed (RH) 


a {S RH } 


O {Y R „} 


Y {T RH } 



Table 2. CEO in scalar- vector-tensor (SVT) classification, i.e., 3- 
tensor notation: 4 x (1+3 + 5) =36 observables. In curly brackets 
we provide an alternative notation of the CEO. 



where we organize the CEO into four groups, which have 
been introduced briefly in the previous section: the active to- 
tal, active handed, reactive total, and reactive handed CEO 
parameter groups, respectively. In addition, we will use a 
coordinate-free 3D formalism and classify the CEO param- 
eters according to rank, i.e., as scalars, 3-vectors, and rank 
two 3-tensors (SVT classification). The 3D CEO are listed 
in Table 12.21 The alternative SVT nomenclature, enclosed 



by cur ly brackets in the the table, was proposed in lOlofsson 
(2001), which used S, V and T, to denote all the scalars, 
vectors, and tensors, respectively. The subscript indexing 
scheme used here is different to that used by Olofsson. The 
CEO 3-tensors are defined as follows. 
The "total" parameters are: 

M = £o r 00 = eo (|E| 2 + |B| 2 ) 12 
P =£ r° = Re [E x B*] /Z 
T =e Q T ij = ul 3 - £o Re [E ig> E : 



-c 2 Bc 



>B* 



(12) 
(13) 
(14) 
(15) 



where I3 is the identity matrix in three dimensions. This 
is the 3D representation of the well-known energy-stress 4- 
tensor e T a P, defined by Eq. ©. 
The "handed" parameters are: 

w =S U 00 = Im [E • B*] /Zq (16) 
V =E U i0 = -Balm [(E x E* + c 2 B x B*)] /2 (17) 
U =e Q U ij = vl 3 - Im [E ® B* - B ® E*] /Z (18) 

This is the 3D representation of the handed energy-stress 4- 
tensor EoU a P, defined by Eq. ©. 
The "reactive total" parameters are: 



/=£ C + =e (|E| 2 -|B| 2 )/2 
-Im[ExB*]/Z 



R =eoQ i0 

1 / 2 
X =e W' Qj0 = - ( e Re [E ® E* - c 2 B ® B*] - -Zl; 



(19) 
(20) 



(21) 



Contrary to the active, total and handed, parameter groups 
above, the reactive total parameter group have no single cor- 
responding 4-tensor. Instead it is composed of parts from 
three different CEO space-time tensors: the vacuum proper- 
Lagrangian defined by Eq. (0, the reactive energy flux den- 
sity from Eq. ( fTol ), and the generalized Stokes parameters 



corresponding to the auto-correlated E and B fields from Eq. 

CD. 

The "reactive handed" parameters are: 
a =£ C_ = - Re [E ■ B*] /Z (22) 
O =£o^4e W = ™Im [(E x E* - c 2 B x B*)] (23) 



Y =eo-e J kl W i0kl = £ ( Re [E $ B* +B ® E*] /Z - ^al 3 



(24) 



Also for this parameter group, there is no single correspond- 
ing 4-tensor. The reactive handed group is composed of parts 
from three CEO space-time tensors: the vacuum pseudo- 
Lagrangian, defined by Eq. ([6j, the reactive handed energy 
flux density from Eq. ([Tol l, and the generalized Stokes pa- 
rameters corresponding to the cross-correlated E and B fields 
from Eq. (fTTT i. 

3 CEO in two dimensions 

Up until now we have assumed that all three Cartesian com- 
ponents of both the electric field, E, and the magnetic field, 
B, are measured. One may ask what happens if some com- 
ponents are not measured; can all the 36 parameters of the 
CEO be retained? Of course this is not possible, some in- 
formation is certainly lost in this case, but what one can do 
is to construct a set of parameters analogous to CEO in two- 
dimensions. As will be shown, a total of 4 x (1 + 1 +2) = 16 
CEO 2D parameters can be derived. 

Assume that we can measure the electric field and the mag- 
netic field in a plane which we can say is the xy-plane without 
loss of generality. Let the two-dimensional (2D) fields in this 
plane be denoted E2D := (E x ,E y ) and B2D := {B x ,B y ), and 
define the scalar product between 2D vectors as 



E2D • B;>d = E X B" X + E y B* 
and the cross product as 

E2D x B 2D = E x B y — E y B x 
and the outer (tensor) product as 



E 2D ®B 



2D 



E X B* E X B* 

EyB* EyB* 



(25) 



(26) 



(27) 



We will however not need to consider all the components 
of the 2D direct product since the 2-tensors we will consider 
are all symmetric and traceless. Hence, we only want the pa- 
rameters which correspond to Pauli spin matrix components 



1 

1 

1 

-1 



(28) 
(29) 



J. E. S. Bergman and T. D. Carozzi: Canonical Electromagnetic Observables 



5 



The Pauli components can be extracted from a 2D matrix by 
matrix multiplying by a Pauli spin matrix and then taking the 
trace (Tr), that is 



Tr{(E 2D ® B* 2 d) <M = (E2D ® B| D ) • • a x 



E X B1 



' EyB* 



(30) 



where the double dots, •■ , have been introduced as a symbo l 
for the double scalar product, see lLebedev and Clo ud ( 2003). 

We can derive a set of two-dimensional canonical electro- 
magnetic parameters from the full CEO by formally taking 

E z =B z = (31) 

and discarding all the parameters that are identically zero. In 
this way we obtain the following set, which we write in the 
2D formalism introduced above. 

The "active total" 2D parameters are: 

-c 2 |B 2D | 2 ) /2 



nil) = £0 ( |E2D 

2 



eo(\E x \*+\E y \- + \cB x 

Re[E 2D xB| D ] /Zo 
Re[E x B;-E y B* x ] /Z 
e Re [E 2D <8>E 



\cB y \ z )/2 



^2D" 
|2 . 



c 2 B 2D ( 



B 



Ta x = 



■■£o(\E x \ 2 -\E y \ 2 + \cB x \ 2 - 
-c 2 B 2D ( 



2DJ 



\cBy\ 



£ Re [E 2D ® E 2D + c z B 2D ® B^ Dj 



o z /2 
/2 

<y x /2 



= eo Re [E x e;+ B x B* y ] 
The "active handed" 2D parameters are: 
w 2 d = Im [E 2D • B 2D ] /Z 
= Im [E x B*+EyB*] /Z 
V z = e Im [E 2D x E 2D + c 2 B 2D x B 2D ] /2 

= e im c 2 b x b;] /2 



U a - 



: Im [E 2D ® B 2D - B 2D <g> E 2D ] ■ • C7 Z /2Z 



= Im [E x B* x -E y B* y ] /Z 
f/^. = Im [E 2D ® B 2D - B 2D ® E 2 rj] • • C7 X /2Z 
= lm[E x B;+E y B x ] /Z 
The "reactive total" 2D parameters are: 

ho 

|2 



£o(|E 2D | 2 -c 2 |B 2D | 2 )/2 



R 7 = 



X, 



eo^r + |^| -\cB x 
Im[E 2D xB^ D ] /Z 
ha[E x B* y -E y B*] /Z 
£oRe[E 2D (8)E5 D 

£o(|^| 2 -|£ v | 2 -|cB,| 2 
£ Re [E 2D ® E 2D - B 2D (gi ■> , 

£ Re c 2 b x b;] 



/2 



-c 2 B 2D ( 



B5o]--ffz/2 

2' 



B 



c»,rj /2 

(7 x /2 



(32) 
(33) 
(34) 
(35) 

(36) 
(37) 
(38) 
(39) 

(40) 
(41) 
(42) 
(43) 



The "reactive handed" 2D parameters are: 

a 2D = Re[E 2D -B| D ] /Z 
= Re[E x B*+E y B* y ] /Zo 

£ Im [E 2D x E 2D - c 2 B 2D x B 2D ] 



O: 



Y a = 



eolm[E x E;-c 2 B x B* y ] 

Re [E 2D ® B 2D + B 2D (» E 2D ] • • (7 Z /2Z 

Ref^B*-^;] /Zo 

Re [E 2D ® B 2D + B 2D (g> E 2D ] • ■ a v /2Z 

Re[£,Z?;+£X] /Zo 



(44) 
(45) 
(46) 
(47) 



We can associate names with these parameters as listed in 
Table [3] The first four parameters, which we call the "to- 
tal" 2D CEO parameters are all well known. These parame- 
ters are also known by different names, e.g., the total energy 
flux is also known as the Poynting vector (z-component), and 
the total energy stress is known as the Maxwell stress tensor 
(difference of diagonal components and off-diagonal com- 
ponent). The remaining three sets of 2D CEO parameters 
are less well known. We will not be able to provide a full 
physical interpretation of each of these parameters; indeed 
their role in space plasma physics is yet to be fully explored. 
We will only mention that the "handed" parameters involve 
spin (helicity, chirality, polarization) weighted energy, i.e., 
the energy of the right-hand wave modes are weighted posi- 
tively and the energy of left-hand wave modes are weighted 
negatively, and these weighted energies are then added. Its 
flux corresponds to the concept of ellipticity and for the case 
of vacuum, it is numerically equivalent to Stokes V param- 
eter. The reactive energy densities come in two groups: the 
"reactive total" and the "reactive handed" 2D CEO param- 
eter groups. From the "reactive total" group, we now rec- 
ognize the reactive energy flux density, as well as the EM 



Symbol 


Name 


"2D 


Total energy 


P z Total energy flux 


Ta z 


Total energy stress a z component 


Ta x 


Total energy stress a A component 


W2D 


Handed energy 


Vz 


Handed energy flux 


Ua- 


Handed energy stress (T.-component 


Uo, 


Handed energy stress C7 A -component 




Vacuum proper-Lagrangian 


Rz 


Reactive energy flux 


X a - 


EM Stokes parameter Q auto-type 




EM Stokes parameter U auto-type 


aiD 


Vacuum pseudo-Lagrangian 


Or 


Reactive handed energy flux 


Ya- 


EM Stokes parameter Q cross-type 




EM Stokes parameter U cross-type 



Table 3. Naming scheme for the 2D CEO parameters. 
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Total epergy flux 



„ 2- 



Tptal stress a r -component 



£ 2- 

1- 



Tptal stress a T -component 



06:45 07:00 07:15 07:30 07:45 
Time in UT (2001-03-31) 



08:00 



06:45 



07:00 



07:15 07:30 07:45 
Time in UT (2001-03-31) 



08:00 




loglO(|xWs/Hz) 
B+2 




■2 
| -4 

1+0.05 



-0.05 
+0.5 


-0.5 

+0.5 


0.5 



□ 





d) 



EM Stokes U param auto-type 
llllUllll.il (l| 



06:45 07:00 07:15 07:30 07:45 
Time in UT (2001-03-31) 



08:00 



Pseudo (vacuum) Lagrangian 



Reactive handed energy flux , 



0.02 



+0.02 


-0.02 



I +0.04 
+0.02 

■-0.02 




Ely) Stokes p param 



pross-type 



; 

fli — [|j 




+0.5 


-0.5 



06:45 07:00 07:15 07:30 07:45 08:00 
Time in UT (2001-033-1) 



Fig. 1. Example dynamic spectra of the 16 normalized two-dimensional CEO parameters. The parameters were computed from STAFF-SA 
data from Cluster space-craft 2 using the ISDAT database system. The following normalization has been applied: each parameter has been 
divided by the total energy except the t otal energy i t self. T hus all spectral values are in dimensionless unit except for the total energy. This 
Figure can be compared with Fig. 1 in lParrot et alj 12003). The 16 parameters are subdivided into a) the active total energy parameters, b) 
the active handed energy parameters, c) the reactive total energy parameters, and d) the reactive handed energy parameters. Note that all the 
quantities are purely electromagnetic in origin and so do not refer to contributions from the plasma. 
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Stokes Q and U parameters, which here are of the auto-type; 
the vacuum proper-Lagrangian needs no further introduction. 
The "reactive handed" group contain the handed counterparts 
of the reactive energy flux density and EM Stokes parame- 
ters, which here are of the cross-typer; the vacuum pseudo- 
Lagrangian is well-known. 



is a physical phenomenon it would be indicative of Faraday 
rotation. Also there seems to be frequency dispersion in the 
handed stress since its components changes sign with fre- 
quency. Finally, the handed energy clearly shows the hand- 
edness of the chorus emissions on its own, without recourse 
to the sign of the total energy flux. 



4 Application of CEO to Cluster data 



5 Conclusions 



Let us demonstrate that the CEO parameters can easily be 
computed from actual data. Assuming that we have mea- 
surements from a vector magnetometer and an electric field 
instrument, all that is required is to auto/cross-correlate all 
measured components and then form the appropriate linear 
combination introduced above. 

As an example we will consid er the STAFF-SA data set on 
the Cluster-II space-craf t missionlEscoubet et al.1 d 19971) . The 
STAFF-SA instrument Icornilleau-Wehrlinl dl997» ~Ts well 
suited for the CEO parameters since it outputs auto/cross- 
correlation of electric and magnetic field components; how- 
ever as Cluster does not measure one of the electric field 
components (namely the component normal to the spin-plane 
of the space-craft) we can only use the 2D version of the CEO 
introduced in the previous section. 

For this particular example, we re-process the high- 
band part of STA FF-SA data from an event discussed in 



Parrot et al.1 d2003l) from 2001-03-31 UT. In Fig. 1 of this 



paper, Parrot et al display certain parameters based on the 
STAFF-SA data computed usi ng a numerical software pack- 
age ca lled PRASSADCO; see lSantolfld(l2003l) : ISantorik et al.l 
(2006). The interesting feature of the 2D CEO parameters is 
that they are the complete set of electromagnetic field ob- 
servables in the spin-plane of the space-craft; and indeed, 
they use up all the parameters in the STAFF-SA dataset ex- 
pect for the magnetic field in the spin direction. Each CEO 
is a distinct physical quantity and examination of the panels 
in Fig. Q] indicates that this is indeed the case, since besides 
showing a common chorus feature (the arch to the left in each 
panel) there are unique points in each of the panels. 

Besides being a complete description of the electromag- 
netic observables, the fact that the CEO parameters are based 
on parameters that conform with the physics of space-time 
means that we can expect physical phenomenon to be mea- 
sured properly. Seeing as how the CEO parameters have 
not been explicitly measured in the past, we can expect that 
their future use may lead to new physical insights, especially 
since several of the parameters are completely new to space- 
physics. As an example consider again the data shown in 
Fig [U It is interesting to note that the reactive total energy 
flux is only significant close to the equator; this implies that 
the equator is the source region for the chorus events, since 
reactive energy flux is typically large close to radiating ob- 
jects due to large standing energy fields. One can also see a 
modulation at 2.5 kHz in the EM Stokes parameters. If this 



The proposed CEO parameters conveniently organize the 
measurements of the full EM wave field. They are physi- 
cally meaningful quantities, i.e. they transform as geomet- 
ric (Minkowski space-time) objects and they are mathemat- 
ically unique (they are irreducible tensors). The CEO retain 
all information, i.e. nothing is lost, and a linear transforma- 
tion back to the full sixtor form exists. Through parameter 
subset selection they could enable considerable data reduc- 
tion. These parameters have clear despinning properties and 
the scalar quantities do not even need despinning. Some of 
the CEO parameters have not been used before to describe 
EM wave fields and can thus be used to reveal new physi- 
cal insights when applied to the analysis of EM wave field 
data measured by spacecraft. To this end a particularly use- 
ful decomposition of the 36 second order EM components 
into twelve 3-tensor quantities have been provided. All the 
CEO are real valued and unique. The active CEO can prop- 
agate to infinity. Notably, the tensors T a P and U a P obey 



the (vacuum) conservation laws d 



respectively, see Bergman et al.l 12008). 



-0 and d a U a P =0, 
The reactive CEO 
on the other hand do not obey conservation laws, and hence, 
can not propagate to infinity. They are nevertheless impor- 
tant since they provide clues to investigate the properties of 
the near field. As illustrated in FigQ] this could be very usful 
for analysing in situ measurements on board satellites when 
they are close to the source. 
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